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Abstract. The performance of a quantum teleportation algorithm implemented on 
an ion trap quantum computer is investigated. First the algorithm is analyzed in terms 
of the teleportation fidelity of six input states evenly distributed over the Bloch sphere. 
Furthermore, a quantum process tomography of the teleportation algorithm is carried 
out which provides almost complete knowledge about the algorithm. 



1. Introduction 

Quantum teleportation [1] is one of the fundamental experiments of quantum 
information science. The transfer of the quantum properties of one system to a second 
(distant) system based on the nonlocal properties of an entangled state highlights 
the most peculiar and fascinating aspects of quantum mechanics. The experimental 
realization of teleportation requires complete experimental control over a system's 
quantum state. For this reason, teleportation has been only implemented in a few 
physical systems [2-9]. One of these systems are strings of cold ions stored in linear 
Paul traps. The achievable level of control over the quantum state of trapped ions 
makes this system an ideal candidate for quantum information processing. Single and 
two-qubit gates constituting the fundamental building blocks for quantum information 
processing have already been demonstrated [10-13] and characterized by quantum 
process tomography [14]. The concatenation of quantum gates in combination with 
measurements has been used for demonstrating simple quantum algorithms [15-17]. 
Quantum teleportation can be viewed as an algorithm that maps one ion's quantum state 
to another ion. In the context of quantum communication, teleportation can also be 
interpreted as a non-trivial implementation of the trivial quantum channel representing 
the identity operation. In this paper, we characterize an ion trap based experimental 
implementation of such a quantum channel by quantum process tomography. We 
improve the previously reported fidelity of the teleportation operation [8] and extend the 
analysis by teleporting the six eigenstates of the Pauli operators a x iV:Z and measuring 
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the resulting density matrices. These data are used for reconstructing the completely 
positive map characterizing the quantum channel. 



2. Teleporting an unknown quantum state 

Teleportation achieves the faithful transfer of the state of a single quantum bit between 
two parties, usually named Alice and Bob, by employing a pair of qubits prepared in a 
Bell state shared between the two parties. The protocol devised by Bennett et al. [1] 
assumes Alice to be in possession of a quantum state ipi n — a|0) + where a and (5 
are unknown to Alice. In addition, she and Bob share a Bell state given by 

|* + > = -L(|0> A |l> B + |l> A |0> B ), (1) 

where the subscripts indicate whether the qubit is located in Alice's or Bob's subsystem. 
The joint three qubit quantum state of Alice's and Bob's subsystem 

\V) AB = -±= (a|00) A |l) B + P\10)a\1)b + a|01) A |0) fl + P\11) a \0)b) (2) 

can be rearranged by expressing the qubits on Alice's side in terms of the Bell states 
= (|10) ± |01>)/V2 and $± = (|00) ± \U))/y/2: 

\*)ab = \(nt<*\l) + 0|O» fl +®a .(«|1> - m)s+n («|0) + P\1))b. 



+y-((3\l)-a\0)) B ). (3) 



By a measurement in the Bell basis, Alice projects Bob's qubit into the states a x ■ 
(o~x °~z) ' ^ini ^in and o z ■ \f in depending on the result of the measurement. If Alice 
passes the measurement result on to Bob, he is able to reconstruct by applying the 
necessary inverse operation of either a x , a z a x , I or a z to his qubit. 

With trapped ions, it is possible to implement teleportation in a completely 
deterministic fashion since both the preparation of the entangled state and the complete 
Bell measurement followed by measurement-dependent unitary transformations are 
deterministic operations. 



3. Experimental setup 

In our experimental setup, quantum information is stored in superpositions of the 
Si/ 2 (m = —1/2) ground state and the metastable D 5 / 2 (m = —1/2) state of 40 Ca + 
ions. The calcium ions are held in a linear Paul trap where they form a linear string 
with an inter-ion distance of about 5 /im. State detection is achieved by illuminating the 
ion string with light at 397 nm resonant with the S1/2 <-> Pi/2-transition and detecting 
the resonance fluorescence of the ions with a CCD camera or a photo multiplier tube. 
Detection of the presence or absence of resonance fluorescence corresponds to the cases 
where an ion has been projected into the \S) or |D)-state, respectively. The ion qubits 
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Figure 1. Teleportation algorithm for three ion-qubits. Ion 1 is prepared in the input 
state \x) = U X \S) while ion 2 and 3 are prepared in a Bell state. The teleportation pulse 
sequence transfers the quantum information to ion 3. During the Bell measurement the 
quantum information in the ions not subjected to the measurement are protected from 
the 397 nm light by shifting the |5)-state population to an additional |D)-substate 
using the pulses denoted by Hide and Hide -1 . The operations labelled X and Z 
represent spin flip and phase flip operations, respectively. 



can be individually manipulated by pulses of a tightly focussed laser beam exciting 
the \S) «-> \D) quadrupole transition at a wavelength of 729 nm. The motion of 
the ions in the harmonic trap potential are described by normal modes, which appear 
as sidebands in the excitation spectrum of the S1/2 -D5/2 transition. For coherent 
manipulation, only the quantum state of the axial center of mass mode at a frequency 
of (jJcom = 27r x 1.2 MHz is relevant. Exciting ions on the corresponding upper or 
blue sideband leads to transitions between the quantum states \S,n) and \D,n + 1), 
where n is the number of phonons. By employing sideband laser cooling the vibrational 
mode is initialized in the ground state |n = 0) and can be precisely controlled by 
subsequent sideband laser pulses. These sideband operations, supplemented by single 
qubit rotations using the carrier transition, enable us to implement an entangling two- 
qubit quantum gate. Further details of the experimental setup can be found in [18]. 

4. Implementing teleportation in an ion trap 

Three ion-qubits are sufficient for the teleportation experiment. One qubit carries the 
unknown quantum information and an entangled pair of qubits provides the necessary 
entangled resource for the information transfer. 

Fig.[Hprovides an overview of the pulse sequence used for teleportation. A complete 
list of all necessary experimental steps is given in Tab. [TJ This pulse sequence can be 
broken down into the following experimental steps: 

(i) Initialization of ion qubits: Initially, the ion string's vibrational motion is 
laser-cooled by Doppler cooling on the S1/2 <-> -P1/2 dipole transition. Subsequent 
sideband cooling on the Si/ 2 «-> -D5/2 quadrupole transition initializes the center-of- 
mass mode in the ground state, which is a crucial prerequisite for the entangling and 
disentangling sideband operations in the teleportation circuit. By a pulse of circular 
polarized 397 nm light, we make sure that all ion qubits are in the Syziyn = —1/2) 
ground state at the beginning of the teleportation sequence. 
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(ii) Bell state preparation: Ion 2 and 3 are prepared in the Bell state (\DS) + 
\SD))/y/2 by a sequence of three laser pulses (see Tab. [I]). We are able to generate 
this entangled state with a fidelity of up to 96% [19]. Furthermore, this particular 
Bell state is highly robust with respect to the major decoherence mechanisms in our 
experimental resulting in a lifetime only limited by the lifetime of the metastable 
D 5 / 2 -level [19]. 

(iii) Preparation of the input state: Ion 1 is prepared in the input state \ipi n ) = 
U X \S), where U x is a single qubit rotation. 

(iv) Rotation into the Bell-basis: In order to carry out the measurement 
in the Bell basis, we have to map the Bell basis onto the product basis 
{|S'S'), \SD), \DS), \DD)}, which is the natural measurement basis in our setup. 
This basis transformation is achieved by first applying a CNOT gate operation to 
the qubits, mapping the Bell states onto separable states, and a final Hadamard-like 
single qubit rotation. In our quantum circuit the CNOT gate, which is extensively 
described in [18], is decomposed into a controlled phase gate and two single qubit 
rotations of length n/2. However, one of the tt/2- rotations (pulse 30 in Tab. [1]) 
is shifted to the reconstruction operations on ion 2. This means that the product 
basis corresponds to a different set of entangled states, namely {($~ + ^ + )/y/2, 
($+ + *-)/V2, ($+ - ^")/v / 2, ($ + - y)/V2} are mapped onto {\DD), \DS), 
\SD), \SS}}. 

(v) Selective read-out of the ion string: Ion 1 and 2 are measured in the product 
basis by illuminating the ions with light at 397 nm for 250 fis and detecting the 
presence or absence of resonance fluorescence on the S\/2 <-» Pi/2-transition that 
indicates whether the individual ion was projected into state \S) or \D). During 
the measurement process the coherence of the target ion 3 has to be preserved. 
Therefore, the S-state population of ion 3 is transferred to an additional Zeeman 
sub-state of the D 5 / 2 level, which is not affected by the detection light [16]. For the 
detection of the fluorescence light of ion 1 and 2 we use a photomultiplier (PMT), 
since its signal can be directly processed by a digital counter electronics which then 
decides which further reconstruction operations are later applied to ion 3. However, 
this requires to read out the two ions subsequently as the states \SD) and \DS) 
cannot be distinguished with the PMT in a simultaneous measurement of both 
ions. This is implemented measuring one ion while hiding the other ion using the 
technique described above. 

(vi) Spin-echo rephasing: Application of the hiding technique to qubit 3 protects the 
quantum information it carries from the influence of the Bell measurement on the 
other ions. However, quantum information stored in the D-state manifold is much 
more susceptible to phase decoherence from magnetic field fluctuations. In order 
to undo these phase errors a spin echo sequence [20] is applied to qubit 3 (pulse 17 
in Tab. [1]). In order to let qubit 3 rephase, a waiting time of 300 fis is inserted 
after completion of the Bell measurement before the reconstruction operations are 
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applied. Simulations of the teleportation algorithm show that the spin-echo waiting 
time which maximizes the teleportation fidelity depends on the chosen input state. 
Since a maximum mean teleportation fidelity is desired, a spin-echo time has to be 
chosen which is the best compromise between the individual fidelities of the input 
states. Additionally, we carry out a spin-echo pulse on ion 1 after the phase gate 
in order to cancel phase shifts during the gate operation. 

(vii) Conditional reconstruction operation: The information gained in step (v) 
allows us to apply the proper reconstruction operations for qubit 3. However, 
compared to the reconstruction operations found in Sec. [2] the preset single qubit 
rotations in our teleportation circuit have to be modified due to the omitted n/2- 
rotation in the Bell measurement and due to the spin echo applied to ion 3 which 
acts as an additional — iF-rotation. First of all an additional 7r/2-rotation is applied 
to ion 3, making up the rotation missing in the Bell analysis. Finally, for the four 
Bell measurement results {\DD), \DS), \SD), \SS)} the single qubit rotations 
{XZ, iX, iZ, 1} have to be applied to qubit 3, i.e. a Z-operation has to be applied 
whenever ion 1 is found in the |D)-state and an X-operation whenever ion 2 is found 
to be in \D). Note that all these single qubit rotations and all following analysis 
pulses are applied with an additional phase <fi. This allows us to take into account 
systematic phase errors of qubit 3, by maximizing the teleportation fidelity for one 
of the input states by adjusting <ft [8]. This optimum phase <fi is then kept fixed 
when teleporting any other quantum states. 

5. Teleportation results 

Due to experimental imperfections and interaction of the qubits with the environment, 
no experimental implementation of teleportation will be perfect. For this reason, we 
describe the experimental teleportation operation by a completely positive map £(p), 
expressed in operator sum representation as [21]: 

4 

£(P) = Yl Xmn A mP A L ( 4 ) 
m,n=l 

where p is the input state to be teleported, and A m e {I, a x , a y , a z } is a set of operators 
forming a basis in the space of single-qubit operators. The process matrix \ contains 
all information about the state- mapping from qubit 1 to qubit 3. 

A useful quantity characterizing the quantum process £ is the average fidelity 
F = J dil)(i/)\£(i/))\i/)) where the average over all pure input states is performed using a 
uniform measure on state space with J dip = 1. In the case of a single qubit process, 
the integral would be over the surface of the Bloch sphere. However, for the calculation 
of F, an average over a suitably chosen finite set of input states suffices [22,23]. Using 
the eigenstates ip±k,k G {x,y,z}, of the Pauli matrices a x ,a y ,a z , F is obtained by 
calculating F = ± £ M±Ei±yi±2} (%|£(^)l%)- 
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Table 1. Sequence of laser pulses and experimental steps to implement 
teleportation. Laser pulses applied to the i-th ion on carrier transitions are 
denoted by Rf(9,(p) and R^(9,ip) and pulses on the blue sideband transition 
by Rf(9, if), where 9 = Qt is the pulse area in terms of the Rabi frequency f2, 
the pulse length t and its phase tp [19]. The index C denotes carrier transitions 
between the two logical eigenstates, while the index H labels transitions from 
the 51/2- to the additional L> 5 / 2 -Zeeman substate used to hide individual ion 
qubits. 





Action 


Comment 




1 


Light at 397 nm 


Doppler preparation 




2 


Light at 729 nm 


Sideband cooling 




3 


Light at 397 nm 


Optical pumping 


angle 


4 


i? 3 + (7r/2,37T/2) 


Entangle ion #3 with motional qubit 


5 


i?C(7r,37r/2) 


Prepare ion #2 for entanglement 


H 


6 


i? 2 + (7T,7r/2) 


Entangle ion 2 with ion 3 




7 

8 


Wait for 1/is - 10 000 /is 

itf(7T,0) 


Stand-by for teleportation 
Hide target ion 


9 




Prepare source ion #1 in state x 




10 


i? 2 + (7r,37r/2) 


Get motional qubit from ion 2 




11 


Rt (7r/v^,7r/2) 


Composite pulse for phasegate 




12 


i?+(7T,0) 


Composite pulse for phasegate 


'55 


13 




Composite pulse for phasegate 


a 


14 




Composite pulse for phasegate 


"33 


15 


itf(7T,7r/2) 


Spin echo on ion 1 


m 

o 


16 


i^(7T,7r) 


Unhide ion 3 for spin echo 


$ 


17 


i^(7T,7r/2) 


Spin echo on ion 3 


,otate 


18 


if (tt,0) 


Hide ion 3 again 


19 


i? 2 + (7T,7r/2) 


Write motional qubit back to ion ^2 


PC 


20 


i?f(7r/2,37T/2) 


Part of rotation into Bell-basis 




21 




Finalize rotation into Bell basis 




22 


i^ 2 (7T, 0) 


Hide ion Z 




23 


r ivi i detection ^pi ^zou /zsj 


Read out ion 1 with photomultiplier 




24 


i??(vr,0) 


Hide ion 1 


+3 

pi 


25 


i? 2 H (7T,7T) 


Unhide ion 2 


Read-< 


26 


PMT detection #2 (250 j*)s 


Read out ion 2 with photomultiplier 


27 


i? 2 H (7r,0) 


Hide ion #2 




28 


Wait 300 /is 


Let system rephase; part of spin echo 




29 


i?3 H (7T,7T) 


Unhide ion 3 


30 


i? 3 c (7r/2,37r/2 + (/)) 


Change basis 


a 

i ° 


31 




ia x "i . conditioned on PMT 


PI +3 
o o 


32 


i? 3 C (7T,7r/2 + 0) 


-\o y ) z detection #1 


Rec 
stri; 


33 


i? 3 c (vr,0) 


ia x conditioned on PMDetection 2 


34 


if (tf x ,^ x + 7T + 0) 


Inverse of preparation of x with offset 


35 


Light at 397 nm 


Read out ion 3 with camera 
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The overlap (i[)j\£(i/;j)\i[)j) between the input state ipj prepared in ion qubit #1 
with the output state generated via teleportation in ion qubit #3 is measured directly 
in our experiment by applying the inverse unitary transformation to ion qubit #3 after 
teleportation and determining the probability to find this qubit in the initial state 15*), 
i.e. formally the teleportation fidelity is given by F te i e = (S\U~ 1 p exp U x \S), where p exp 
is the quantum state of ion qubit #3 after teleportation. For the six input states 
^ = \S), V> 2 = \D), ^3 = (\D) - i\S))/y/2, ^4 = (\D) - \S))/y/2, ^ = (\D)+i\S))/y/2, 
ipe = (\D) + \S))/y/2, the teleportation fidelities range between 79% and 87% (see Fig. 
12]), with an average fidelity of F — 83(1)%. This average fidelity proves successful 
operation of the teleportation algorithm, as it exceeds the maximum value of 2/3 that 
is achievable without using entangled states [24]. 
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Figure 2. Measured teleportation fidelity -Fi e ; e for six different input states. All 
fidelities are well above the 2/3 threshold proving successful quantum teleportation. 
The average teleportation fidelity is F te i e = 83(1)%. 



A more complete way of characterizing the teleportation process is achieved by 
determining the output state of qubit 3 by quantum state tomography, which requires 
measurements in three different measurement bases. From these measurements, the 
density matrix of the output qubit is estimated using a maximum likelihood algorithm 
[19]. The resulting density matrices of the six input states are shown in Fig. [3] 

Full information about the relation between the input and output of the 
teleportation algorithm is gained by a quantum process tomography. This procedure 
requires to determine the output state £(pi) after application of the investigated 
operation for a set of at least four linear independent input states pj. With this data, 
the process matrix x is obtained by inverting equation (j4]). Due to inevitable statistical 
errors in the measurement process the resulting x will m general not be completely 
positive. This problem is avoided by employing a maximum likelihood algorithm, which 
determines the completely positive map which yields the highest probability of producing 
the measured data set. We use the tomographically reconstructed input states ijJx-ijjQ 
for a determination of the process matrix x by maximum likelihood estimation [25]. 
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Figure 3. Real and imaginary part of the density matrix of the output qubit for 
the six different input states a) xpi = \S), b) ip 2 = \D), c) ip 3 = (l-D) - i\S))/^2, d) 
^ 4 = (\D) - \S))/V2, e) = (\D) + i\S))/V2, f) A = (\D) + \S))/y/2. 



The absolute value of the elements of the resulting process matrix Xteie is shown in Fig. 
Hk). As expected, the dominant element is the identity with xn — 0.73(1), which is 
identical to the process fidelity F proc = tr(xidteieXteie), where Xidteie denotes the ideal 
process matrix of the teleportation algorithm. This agrees well with the average fidelity 
stated above, as average and process fidelity are related by F = {2F proc + l)/3 for a 
single qubit map [23]. 

A quantum process operating on a single quantum bit can be conveniently 
represented geometrically by picturing the deformation of a Bloch sphere subjected 
to the quantum process [21]. The quantum operation maps the Bloch sphere into itself 
by deforming it into an ellipsoid that may be rotated and displaced with respect to 
the original sphere representing the input states. This transformation is described by 
an affine map r out = OSr in + b between input and output Bloch vectors where the 
matrices O and S are orthogonal and positive-semidefinite, respectively. Fig. Hb) shows 
the result for the teleportation algorithm. The transformed ellipsoid is centered at 
b ~ (0, 0.09, —0.05) with errors of about ±0.03 for each coordinate. The matrix S 
shrinks the sphere anisotropically (its eigenvalues are 0.78, 0.58, 0.55), O rotates the 
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Figure 4. Results of process tomography of teleportation algorithm. In (a) the 
absolute value of process matrix x is shown. The dominating diagonal element is the 
identity with \ii — 0.73(1). The plot drawn in (b) shows how the input states lying 
on the surface of the initial Bloch sphere (meshed surface) are transformed by the 
teleportation algorithm, with the output states lying on the solid surface. 

sphere by an angle of about 2° (2). This demonstrates that the loss of fidelity is mostly 
due to decoherence and not caused by an undesired unitary operation rotating the sphere 
as the orientation of the deformed Bloch sphere hardly differs from the orientation of the 
initial sphere. The results are consistent with the assumption that the rotation matrix 
O is equal to the identity as desired. 

6. Conclusion 

We demonstrated deterministic teleportation of quantum information between two 
atomic qubits. We improve the mean teleportation fidelity F = 75% reported 
in [8] to F — 83(1)% and unambiguously demonstrate the quantum nature of the 
teleportation operation by teleporting an unbiased set of six basis states [26] and 
using the data for completely characterizing the teleportation operation by quantum 
process tomography. The process tomography result shows that the main source of 
infidelity is decoherence while systematic errors are negligible. To make further progress 
towards high-fidelity quantum operations, decoherence rates have to be reduced by 
either reducing environmental noise or encoding quantum information in noise-tolerant 
quantum states [27]. 
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